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We investigate high-Reynolds number turbulence in dilute polymer solutions. We show the exis- 
tence of a critical value of the Reynolds number which separates two different regimes. In the first 
regime, below the transition, the influence of the polymer molecules on the flow is negligible and 
they can be regarded as passively embedded in the flow. This case admits a detailed investigation of 
the statistics of the polymer elongations. The second state is realized when the Reynolds number is 
larger than the critical value. This regime is characterized by the strong back reaction of polymers 
on the flow. We establish some properties of the statistics of the stress and velocity in this regime 
and discuss its relation to the drag reduction phenomenon. 

PACS numbers 83.50.Ws, 61.25.Hq, 47.27.-i, 05.40.-a 
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I. INTRODUCTION 

In this paper we present a theoretical investigation 
of turbulence in dilute polymer solutions. As opposed 
to Newtonian fluids, such solutions possess additional 
macroscopic degrees of freedom related to the elasticity 
of the polymer molecules. Relaxation times of elastic 
stresses can be comparable with time-scales of the flow 
which means that the relation between the stress and 
velocity gradient is non-local. Due to the non-trivial in- 
teraction of inertial and elastic degrees of freedom, the 
polymer solutions exhibit a variety of regimes. For ex- 
ample, a turbulent-like state has been recently observed 
at very low Reynolds' numbers 0]. Here we will con- 
sider the more familiar situation of turbulence at high 
Reynolds numbers. Probably, the most striking effect of 
polymers on the high Reynolds number flows is the drag 
reduction phenomenon. The addition of long-chain poly- 
mers in concentrations as small as 10~^ by weight can 
induce a substantial reduction of the drag force needed 
to push a turbulent fluid through a pipe [^|-|| . 

The reason why small amounts of polymer can signif- 
icantly modify properties of the fluid is the flexibility of 
polymer molecules. At equilibrium a polymer molecule 
coils up into a spongy ball of a radius Rq. The value of 
i?o depends on the number of monomers in the molecule, 
which is usually very large. For a dilute solution with 
the concentration, n, satisfying uRq ^ 1, the influence 
of equilibrium size molecules on the hydrodynamic prop- 
erties of the fluid can be neglected. When placed in a 
non-homogeneous flow, such a molecule is deformed into 
an elongated structure that can be characterized by its 
end-to-end distance R. If the number of monomers in a 
typical polymer molecule is large, the elongation R can 
be much larger than Rq. The influence of the molecules 
on the flow increases with their elongation and may be- 
come substantial when R^ Rq. 



The deformation of the molecule is determined by two 
processes, the stretching by the velocity gradients and 
the relaxation due to the elasticity of the molecule. To 
understand how a molecule resists the deformation by 
the flow, let us consider its relaxation. Recent experi- 
ments with DNA molecules indicate that the relaxation 
is linear in the wide region of scales Rq <^ R <^ i?max, 
where i?max is the maximum extension . In the case of 
polymer molecules, theoretical arguments and numerics 
presented in support the linear relaxation. These re- 
sults can be understood if we assume that at i? 3> -Rq the 
role of excluded volume and hydrodynamic interactions 
between the monomers becomes negligible. Then the ran- 
dom walk arguments suggest that the entropy of polymer 
molecules is quadratic in the range Rq 'ti R ^ Rm&x im- 
plying linear relaxation. Whether the polymers are ex- 
cited by the flow is determined by the softest relaxation 
mode that describes the dynamics of the elongation R. In 
the absence of stretching, the relaxation of R is described 
by R — ~R/t, where t is the largest of the polymer re- 
laxation times. The relaxation time is i?-independent at 
Rq ^ R ^ ^max- If the end-to-end distance R is of the 
order of the maximum extension, r starts to depend on 
R and the dynamics of the molecule becomes nonlinear. 

The behavior of the molecule in an inhomogeneous 
steady flow depends on the value of the Weissenberg num- 
ber, Wi, deflned as the product of the characteristic ve- 
locity gradient and r. When a polymer molecule is placed 
in a flow, smooth at the scale R, the velocity difference 
between the end-points is proportional to R multiplied by 
the characteristic value of velocity gradient. At Wi ^ 1 
relaxation is fast compared to the stretching time and the 
polymer always relaxes to the equilibrium size, Rq. The 
behavior of the polymer at Wi ^ 1 depends on the geome- 
try of the flow. For purely elongational flows the molecule 
gets aligned along the principal stretching direction. If 
the velocity gradient is larger than the inverse relaxation 
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time, i.e. Wi > 1, the elastic response becomes too slow 
in comparison with the stretching and the molecule gets 
substantially elongated . The sharp transition from the 
coiled state to the strongly extended state is called the 
coil-stretch transition. Rotation can suppress the transi- 
tion and even damp it completely since the molecule does 
not always point in the stretching direction (see e.g. [^). 
For example, no coil-stretch transition occurs in the case 
of a shear flow, which is a combination of elongational 
and rotational flows. 

In contrast to the steady flows, a polymer molecule 
moving in a smooth random flow alternately enters re- 
gions of high and low stretching. As the intensity of the 
flow increases the effect of the stretching becomes more 
pronounced. One can generally assert the existence of the 
coil-stretch transition. This has been flrst demonstrated 
by Lumley who considered the situation where the 
characteristic time of variations of velocity gradient is 
much larger than the inverse of the characteristic value 
of the gradient. He showed that if the amplitude of ve- 
locity gradient fluctuations is large enough, the expec- 
tation value of grows with time, which signifies the 
coil-stretch transition. We have demonstrated in |^ that 
the coil-stretch transition occurs in any random flow and 
established a general criterion for the transition. In par- 
ticular, the transition occurs in the situation where the 
time of velocity gradient variation is of the order of the 
inverse of its characteristic value, which is likely to be 
the case for real flows. The coil-stretch transition in ran- 
dom flows is controlled by the parameter Air, where Ai 
is the average logarithmic divergence rate of nearby La- 
grangian trajectories, to be referred to as the principal 
Lyapunov exponent. It is positive for an incompress- 
ible flow [ p^|jTl| ]. The molecules are weakly stretched if 
AiT < 1 and strongly stretched otherwise. Therefore for 
random flows the parameter Air plays the role of the 
Weissenberg number. 

To describe the behavior of a polymer molecule in tur- 
bulent flows, let us briefly review the basic properties of 
turbulence of incompressible Newtonian fluids. A high 
Reynolds number flow consists of chaotic motions from 
a wide interval of scales, rj <^ r ^ L, where L is the 
scale at which the flow is excited and rj is the viscous 
scale. The energy pumped at the scale L cascades down 
to the scale rj, where it is dissipated. The size of polymer 
molecules is usually much smaller than the viscous scale. 
Viscosity makes the flow smooth at r <C ?7, i.e. the veloc- 
ity difference between two points is given by the velocity 
gradient multiplied by the distance. Then, the stretch- 
ing of molecules is determined by the gradient of veloc- 
ity, which should be considered random in a turbulent 
flow. The Lyapunov exponent can be estimated as the 
characteristic value of the velocity gradient, which is de- 
termined by the eddies at the viscous scale. As Reynolds 
number increases, the velocity gradient increases, and so 
does Air. At some value of Reynolds number it reaches 
the value 1 and the coil-stretch transition occurs. 

Several mechanisms can limit the stretching of poly- 



mers. The first one is the internal non-linearity of the 
elasticity of the polymer molecules. If this mechanism 
dominates, then above the transition the molecules are 
stretched up to the maximal elongation -Rmax- An al- 
ternative explanation has been proposed by Tabor and 
de Gennes It is based on the fact that if the elon- 
gation of a polymer molecule is larger than the viscous 
length of turbulence, rj, the elastic force always wins over 
the stretching. Estimates using the parameters of typ- 
ical polymer solutions show that this situation is diffi- 
cult to realize. Therefore, we will assume that the in- 
equality i? ^ 77 is satisfied. It will also enable us to 
write local equations describing the dynamics of elastic 
stresses. Another mechanism is the back reaction of the 
polymers on the flow. It is caused by the collective con- 
tribution of coherently deformed polymer molecules into 
the stress tensor. This elastic part of the stress grows 
with the elongation of the molecules. When it becomes 
of the order of the viscous stresses existing in the flow, the 
polymers modify the flow around them and the stretch- 
ing diminishes. As a result, a dynamic equilibrium is 
realized at the characteristic elongation, -Rback- The to- 
tal polymer stress is proportional to nR^, so that -Rback 
depends on the polymer concentration n. Therefore if 
the concentration is large enough, the value of i?back is 
much smaller than i?max- We will consider the effect of 
the back reaction under the assumption i?back <C -Rmax- 
Probably, the condition i?back ^ -Rmax is necessary for 
the existence of a stationary state, because the polymer 
molecules stretched up to -Rmax are intensively destroyed 
by the flow. 

Above the coil-stretch transition the back reaction 
modifies the small-scale properties of turbulent flows, 
which leads to the emergence of a new scale, > rj. 
Large-scale eddies with the sizes r > do not excite 
elastic degrees of freedom, so that the usual inertial en- 
ergy cascade is realized at these scales. At smaller scales 
inertial and elastic degrees of freedom exchange energy, 
which is dissipated mainly due to the polymer relaxation. 
The energy cascade terminates at , so that plays the 
role of a new dissipation scale. 

The plan of the paper is as follows. In Sec. || we 
introduce a system of equations describing the coupled 
dynamics of inertial and elastic degrees of freedom. In 
Sec. Ill we study the situation when the back reaction 



of polymers is small and can be disregarded. We find the 
probability distribution function of the elastic stress ten- 
sor and examine its correlation functions. In Sec. IV we 
study the influence of the back reaction on the flow and 
establish some properties of the velocity and stress statis- 
tics. In Conclusion we summarize our results and discuss 
their implications for the drag reduction phenomenon. 
In appendix ^ we present a detailed derivation of the 



probability density function from Sec. Ill, Appendix 



is devoted to a simple model illustrating some aspects of 
the interaction between the flow and polymers. Prelimi- 
nary results of this work have been published in M . 
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II. BASIC RELATIONS 

Following Hinch let us consider the dynamics of a 
polymer molecule in a smooth velocity field. The degree 
of freedom related to the elongation of the molecule is 
described by the vector R, connecting the end-points of 
the molecule. The equation describing the dynamics of 
R in the absence of a surrounding flow is 



^tR^ 



(2.1) 



where E is the free energy of the molecule, Q is the ther- 
mal noise, and T is the kinetic coefficient which deter- 
mines the relaxation of the molecule. The correlation 
function of C is 



(2.2) 



where T is temperature and ks is the Boltzmann con- 
stant. If the size of molecules is much smaller than the 
viscous length, which we assume, the molecule moves in 
a constant gradient flow. Its influence is described by the 
equation 



^tR^ - RjVjv, + r— - , 

Otii 



(2.3) 



where the stretching term —RjVjVi is added to Eq. (2.1) 



The velocity derivative must be evaluated at the position 
of the molecule. To avoid misunderstanding note that we 
mean the flow which is "external" to the molecule, ex- 
cluding the velocity induced by the relative motions of 
its chains. 

The entropy of the molecule has a quadratic depen- 
dence on the elongation, R, in a wide interval 0. It 
implies that the molecule can be treated in terms of elas- 
ticity theory with a Hook modulus Kq, so that E — 
KqR^ /2. This expression is correct provided R ^ -Rmax, 
where -Rmax is the maximum elongation of the molecule. 
The equilibrium size of the molecule, Rq, can be esti- 
mated from the condition E ^ ksT as yJksT / K^. Sub- 
stituting the energy into Eq. (p. 3D wc get 



dtR,. 



+ RJt = 



(2.4) 



where r = {TKq)"^. We see that r is the molecular 
relaxation time. 

Generally, the kinetic coefficient F or the relaxation 
time r is a function of R, which reflects non-linear charac- 
ter of the molecule relaxation ||l5| , related to such effects 
as internal hydrodynamic interaction of chains in the 
polymer molecule. For example, the flnitely extendible 
nonlinear elastic model |15| assumes r cx 1 — R? / R^^^. 
One expects that a dependence of F and t on R can be 
disregarded at i? <C -Rmax- Below we study this situation. 
Possible statistical consequences of the non-Hookean de- 
pendence of the free energy on R have been investigated 
in ra. 



Equations (2.3,2.4) assume that R is the only mode 
related to the molecular deformation, which is an ide- 
alization. Actually, the molecule has many deforma- 
tional degrees of freedom that have different relaxation 
times. They have been observed experimentally Q . Nev- 
ertheless, in the turbulent flows, only the mode with the 
largest relaxation time is strongly excited whereas other 
modes are excited at most weakly. Thus, Eq. (2^) should 
be considered as the equation describing the principal 
mode. 



A. Continuous media equations 

To study the dynamics at the scales much larger than 
the inter-polymer distance the polymer solution can be 
regarded as a continuous medium. The appropriate de- 
scription is done in terms of macroscopic quantities which 
are averages of microscopic variables over the volume. 
The polymer molecules are characterized by the average 
conformation tensor 



A, 



(2.5) 



The volume of averaging should contain a large number 
of polymer molecules and be smaller than the charac- 
teristic scales of the processes under consideration. The 
tensor A can also be interpreted as the average over the 
statistics of the the rmal noise C- Equation for Aj^ follows 
fromEqs. (p^,|2^) 



9tAy + {vV)A,j = AkJ\'kV^ 
2 

— (Aij - AoSij) , 



Aik'^kVj 



(2.6) 



where Ag — ksT/KQ. Equation (2.6) is linear in A, 
which is correct provided A <^ Ramx- 

Equation should be supplemented by the equa- 

tion for the fluid velocity. This equation is a consequence 
of the momentum conservation law. In order to derive it 
one should take into account the contribution of the inner 
elastic forces of polymer molecules to the total stress ten- 
sor of the fluid. If n^A; is the elastic stress tensor per unit 
mass, the polymer contribution is gilik where g is the 
mass density of the fluid. In the Hookean approximation 



Kpn 
Q 



Aik - UqS; 



ik ' 



(2.7) 



where Hq = KoAon/ g = ( n/ g)kBT originates from the 
thermal noise ^ in Eq. (2.4). Here n is the concentration 
of the polymer molecules. If the flow is incompressible, 
i.e. p is a constant and \7 ■ v — 0, the momentum conser- 
vation law reads 

dtv, + {vV)v, + g'^V.P = vV^v, + VfeH.fc + /, . (2.8) 

Here P is the pressure, v is the kinematic viscosity of 
the solvent and / is the e xtern al force (per unit mass) 
driving the flow. Equation (2.8) is a generalization of the 



3 



Navier-Stokes equation to the case of viscoelastic fluids. 
To simplify the consideration we assume that / is ho- 
mogeneously distributed over space. It is a common be- 
lief that this case does not differ qualitatively from that 
realized for real experimental setups, where pumping is 
usually related to the boundaries. 



The applicability condition of Eq. (2.6) is A <C -Rmax- 
It can be rewritten in terms of elastic stress tensor as 

n < n 

max — I^qtiq "'^-R^ax' assumc that this condi- 
tion is satisfied for relevant fluctuations. The interaction 
between v and 11 turns on if 11 exceeds the viscous stress, 
i/Wv. The latter can be estimated as i^Xi where Ai is the 
average logarithmic divergence rate of nearby Lagrangian 
trajectories. Under the condition Umax ^ i^^i the poly- 
mer molecules exert no influence on the flow except for 
a small renormalization of the viscosity of the solution. 
Thus, the inequality Umax 2> vXi is a necessary condition 
for the polymers to have non-trivial effects. 

The free energy of the viscoelastic fluid is the sum of 
the kinetic and elastic contributions 



dr 



8 2 

-V 

2 



[tr A - Aq ln(det A/ Aq 



where the second term represents the entropy of the 
molecules. Then we find from Eqs. (2.6-2. 



'dt 



g dr f v 



Q dr I vi^iVjY + -tr U (Hq + U) U 



(2.9) 



This equation provides a mathematical formulation of 
the energy balance: the force supplies the energy, which 
is then dissipated due to viscosity and relaxati on o f poly- 
mers. Note that the second integral in Eq. (2^) has a 
definite sign, as it should be for a dissipative term. Rela- 
tive contributions of the viscous and the elastic terms to 
the energy dissipation can be different. 

If the forcing is statistically homogeneous, then statis- 
tically homogeneous steady state is realized. It can be 
described in terms of correlation functions of v and A 
(or V and 11), which are averages over the statistics of 
the pumping force / or over space. In the steady state 
the average value of dT jdt is equal to zero. Therefore 
we get from Eq. (2.9) 
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v(X^,v,f) + - (tr [n (n + no)-^n] ) = e , 



(2.10) 



where e = (/'u) is the mean energy injection rate (per 
unit mass) by the external force. 

Generally, the diffusion term kV j4_ should also be 
added to the right-hand side of Eq. (2.6). The diffusion 
coefficient, k, is small due to a large number of monomers. 
It is possible to show that the limit k ^ is regular, so 
we can disregard the diffusion. The diffusion term can 



play a minor role for scales r < ^ k/\i. If \J k/Xi is 
smaller than the intermolecular distance, then the diffu- 
sivity is irrelevant in the whole region of applicability of 
the macroscopic approach. 



B. Lagrangian Description 



Equation (2.6) can formally be solved in the La- 
grangian reference frame. Let us introduce A{t, r) = 
A[t,x{t,r)], where x{t,r) is the Lagrangian trajectory 
defined by the relations 



dtx = v{t,x), x{ta,r) = r. 



(2.11) 



The condition x{to,r) =r ensures that the fields A and 
A coincide a.t t ^ t^. The point r plays the role of a La- 
grangian marker. The tensor A(t , r) satisfies the matrix 
equation following from Eq. ( |2.6| ) 

atA = aA + Acr^- ^ (a- Ao) , (2.12) 
a.,j{t,r)^VjV.,[t,x{t,r)]. (2.13) 

Here a is the tensor of the velocity derivatives in the 
Lagrangian frame and the superscript T denotes a trans- 
posed matrix. Due to the causality, the value of the field 
A{t, r) is determined by its dynamics at times t' < t. 
Therefore we will be interested in the backward in time 



Lagrangian evolutio n of A described by Eqs. (2.11,2.12). 

A solution of Eq. ( 2.12 ) can be written in terms of the 
Lagrangian mapping matrix W defined by the relations 



dtW{t,t') = c7{t)W{t,t') , W{t',t') = l. (2.14) 

The matrix W is defined for a given Lagrangian tra- 
jectory, and therefore it depends on its marker r. For 
brevity we omit the argument r in W. The matrix W 
describes the deformation of infinitesimal fluid volumes. 
For example, the separation, 6x, between two close La- 
grangian particles changes according to 



5x{t) = Wit,t')6x{t') . 



(2.15) 



( pT5| ) that W^J{t,to,r) 



It follows from Eq. 

dxi{t,r)/drj. The incompressibihty condition V • t; = 
is formulated in terms of cr as tr cr = 0. A consequence 



of Eq. (2.14) is 



det = 1 . 
Using Eqs. ( |2.12| , ^.14D we obtain 
2Aoe / 



(2.16) 



A{t,r) 



dt'W{t,t',r)W^{t,t\r)e^''/' 



At t = to this equation gives 
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A{to, r) = ^ / dt W{to,to - t, r)W^{to, k - t, r)?,-^*/- . 
T Jo 



(2.17) 



It is easy to understand the meaning of Eq. ( 2.17 ). The 
polymers are advected along the Lagrangian trajectories 
being stretched by the velocity gradient and relaxing to 
their equilibrium shape due to elasticity. The value of the 
conformation tensor A is determined by the sum of the 
contributions of these processes at earlier times picked 
along the Lagrangian trajectory arriving at r. The term 
WW'^ describes the stretching and the exponential fac- 
tor accounts for relaxation. 



Expression ( 2.17 ) shows that when calculating corre- 
lation functions, the volume averages can be substituted 
by averages over the statistics of W. 



III. PASSIVE ELASTIC TENSOR IN A RANDOM 
VELOCITY FIELD 

In this section we consider the polymer molecules as 
passive objects, i.e. we assume that the inertial proper- 
ties of the fluid are not perturbed by the polymer elas- 
ticity. In ot her words, we assume that the term Vfcllife 
in Eq. (2.S) can be disregarded so that the dynamics 
of V is independent of the polymer dynamics. Equation 
( ^.6|) determines the statistics of the conformation ten- 



sor A provided the statistical properties of v are known. 
We consider a high Reynolds number flow and assume 
that its statistics is stationary, spatially homogeneous 
and isotropic. 

One might think that in order to examine the correla- 
tion functions of A one needs to know the precise statis- 
tics of the velocity field described by the Navier-Stokes 
equation. However, we will show that the statistics of A 
is universal, i.e. it does not depend on the details of the 
velocity statistics. The crucial property underlying the 
universality is a finite Lagrangian correlation time of the 
velocity deri vati ves matrix a Jl7[ . 

Equation ( |2.6| ) shows that the Lagrangian dynamics of 
the polymer stress tensor is determined by velocity gra- 
dients, which are related to the viscous scale r] of the tur- 
bulence. A typical value of the velocity gradient can be 
estimated as Ai, where Ai is the logarithmic divergence 
rate of nearby Lagrangian trajectories. Then it follows 
from Eq. (2.8) that the feedback of the polymers on the 
flow can be neglected if 11 ^ lyXi. This is the applicabil- 
ity condition of the passive approach. As we show below, 
the condition 11 <^ vXi is satisfied for typical fluctuations 
if AiT < 1. 



A formal solution of Eq. (2^) is given by Eq. ( 2.17 ). 
The condition Ait < 1 means that the exponentially de- 
caying factor exp(— 2i/r) in Eq. ( ^.17 ) dominates over 
the product WW^ ~ exp[2Aii]. In this case, for a typ- 
ical velocity fluctuation the integral over t converges at 
t ~ r and therefore A fluctuates near Aq. In addition to 
the strong peak at A ^ Aq, the probability distribution 



function of A has a power tail. This has been demon- 
strated in 1^ in terms of the molecule elongation R. The 
fluctuations A Aq are formed if the product WW^ in 
Eq. ( 2.17 ) is anomalously large for a long time. Such 



events can be described in terms of a universal scheme 
(see Ref. [Esl) shortly presented below. 



A. Long-Time Lagrangian Statistics 

Let us briefly review the long-time statistical proper- 
ties of th e Lagrang ian mapping matrix W, determined 



by Eqs. ( ^.13|j2.14| ). We consider W{t+,t_) at t+ > t_ 
and assume that — i_ is much larger than the La- 
grangia n cor relation time To- of the velocity derivatives 
matrix ( ^.13| ) . If the velocity statistics is homogeneous in 
time, the probability distribution of W(t+ A^) depends 



on the difference i+ — t_ only. Equation ( 2.14 ) implies 
that at t+ — ^ the matrix W \s & product of a 
large number of independent matrices. This is the main 
reason for the universality in the statistics of W . 
It is convenient to decompose the matrix W as 



W{t+,t^) = MKN , 



(3.1) 



where A is a diagonal matrix, and M and N are orthog- 
onal matrices ll9j. We denote the diagonal elements of 
A as e^^ , e''^ , and e^^ , and assume that they are or dered : 
pi > P2 > P3- As a consequence of t he co nstraint (2.1i) 
we have pi+ P2+ Pi — 0. Equation (^.14 ) can be rewrit- 
ten in terms of p^, and the matrices M and N. The 
equations for pi are 



dpi 
dt+ 



(3.2) 



where a = M'^aM and no summation over the repeat- 
ing index i is implied. The matrices M and N satisfy 
dtN = riiN and dtM = A/f^a, where 



2sinh(pi - Pk) ' 



^2ik — 



CTike 



2pk 



2pi 



for i ^ k and = ^2ik = for i — k. It is possi- 
ble to show that the eigenvalues of W repel each other, 
so that the inequalities e^^ ^ e''^ ^ e^^ are satisfied at 
t+ — t- ^ To- |lj]. Then the matrix fli tends to zero 
exponentially fast, i.e. N is determined by times of the 
order of in the vicinity of t_. The matrix 1^2 becomes 
/9-independent at t+ — t_ ^ To- and the evolution of M 
is decoupled from that of pi. Then the value of M is 
determined by the time of the order of at t w i.e. 



at ^4 



t- > 



The solution of Eq. (3.2 



it becomes t_ 
is 



■independent. 
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(3.3) 



where the right-hand side of Eq. (3.3) is an integral of a 
random process independent of pi. Equation (|3.3|) shows 
that the variables pi fluctuate around their average val- 
ues Xi{t+ — t-). Here the constants are equal to {da). 
They are called the Lyapunov exponents of the flow. 
Generally, the spectrum of the Lyapunov exponents is 
non-degenerate: Ai > A2 > A3, which is a necessary con- 
dition for the formalism to be self-consistent. The incom- 
pressibility condition ensures the identity Ai-|-A2-I-A 3 — , 
which implies Ai > and A3 < 0. Using relation ( 2.15 ) 
one can show that Ai is indeed the average logarithmic 
rate of the divergence of two Lagrangian trajectories: 



d\n\5p\\^ 



dt 



Ai 



Similarly, Ai -I- A2 = — A3 is the average logarithmic rate 
of the area growth. 

Under the condition t+ — t_ 3> tv the quantity pi can 
be considered as a sum of a large number of indepen- 
dent random variables. It is known from the statistical 
mechanics (see e.g., |^) that the distribution of such 
quantities is given by the exponent of an extensive func- 
tion. In our case the probability distribution function 
(PDF) of p, is 



'P{t,Pi,P2,P3) oc iexp 
xd{pi + P2 + ps) , 



-tS 



Pi - Alt p3 - X^t 



t 



t 



(3.4) 



where t ~ — and pi > p2 > Pz is implied | |i8| |. The 
main exponential factor of the PDF has a self-similar 
form described by the function S called entropy function 
(see |l8| , ^l| , p2[ ) . It is positive, convex and has a minimum 
at zero values of its arguments. The precise form of S 
is determined by details of the velocity statistics. The 
PDF has a sharp maximum at pi — Xit. In its vicinity 
the function S has a quadratic expansion, i.e. the dis- 
tribution of p is Gaussian. However, if one is interested 
in the expectation values of exponential functions of pi, 
they are determined by the wings of the PDF where the 
Gaussian approximation is invalid. This entails the use 
of the whole entropy function. 

To average the functions of pi only, one can introduce 
the reduced probability distribution function. 



V{t,pi) (X ■^f'xp 



-tSi 



pi - Alt 
t 



(3.5) 



which is an integral of V{t, pi, P2, Pz) over p2 and p^. At 
small X the function Si{x) can be written as 



(3.6) 



Here A = j dt {{(Tii{t)dii{Q))) (where double brackets 
designate irreducible correlation function) determines the 



dispersion of pi: {{pi — XitY) ~ tA. Expansion (3^) is 
sufficient to describe typical fluctuations of pi, whereas 
the whole function S is needed to describe rare events. 

In the passive regime the statistics of velocity gradients 
is determined by the fluctuations at the viscous scale rj. 
The Lagrangian correlation time is the turnover time 
at this scale. It can also be estimated as A^^. Using the 
expression e = iy{{Wv)'^) for the energy dissipation rate 
one can write the estimates Ai ~ \f^Jv and A Ai for 
the Lyapunov exponent and the dispersion. 



B. Single-Point Statistics 

In this subsection we examine the single-point statis- 
tics of the conformation tensor A at Ait < 1. As ex- 
plained above, most of the time A fluctuates near Aq. 
We are interested in large values A Aq because it is 
only for large values of A that the polymers can possibly 
lead to noticeable effects. Large values are determined 
by the velocity fluctuations such that the product WW^ 
is anomalously large for a long time. To find t he t ail of 
the PDF of A let us substitute decomposition (|3.lD into 
Eq. (12.171). We obtain 



M^AM = ^ dtA\t)c^p 
r Jo 



2t 

T 



where we used the t-independence of M at large t. Under 
the condition e^'^ ^ e^^ ^ e^^ the tensor A is uniaxial: 



A^j ~ TniTij 



(3.7) 



Here n is a unit vector, ni ~ Mn, uniformly distributed 
over a sphere, and T = tr A: 



2Ao f°°,^ 
/ dt exp 

T Jo 



2pi(t)-- 

T 



(3.8) 



Thus the matrix A is expressed via the scalar T, which is 
independent of p2 and p^. The statistics of T cannot be 
directly examined in terms of the single-tim e pr obability 
distribution function (3^) because integral ( ^ ) involves 
diffe rent times. Nevertheless, it is possible to use PDF 
( |3.5[ ) to find the asymptotic behavior of the PDF of T at 
T 3> Aq. a rigorous derivation is presented in Appendix 
(cf. [HH]). Below we use a simple semi-qualitative 
method. 



For a typical fluctuation of the velocity, integral (3.S) 
converges at t ~ r, which gives T ^ Aq. To find the prob- 
ability distribution for large deviations of T one should 
analyze rare events leading to a given value T ^ Aq, 
and find the event with the maximum probability. Let 
us establish the structure of such fluctuations. It is obvi- 
ous that pi should initially grow faster than t/r during 
some interval of time. To ensure non-zero value of the 
probability of such configuration, pi should then return 
to its average value, Xit. Since Ai < 1/t, the difference 
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Pi{t) — t/r has a maximum at some time, t = t^,. At 
T ^ Ao the maximum is sharp and integral (3_^) is de- 



termined by its vicinity. With the logarithmic accuracy 



ln(T/Ao)«2pi(t,)-2t,/r. 



(3.9) 



The probability of the event is also determined by a vicin- 
ity of pi(i*) because it corresponds to the maximal de- 
viati on f rom the average value of pi . In accordance with 
Eq. (3^) it can be estimated as 



InP « -t^Si 



-Ai 



Substituting here pi{t^) expressed via T from Eq. (3.£) 
and maximizing the result over we get the condition 



5i Ai^ -/35; (^/3+i-Ai^ =0, (3.10) 

where (3 = (2t*)^^ ln(T /Ao) . Using the convexity of 5*1, 
one can show that Eq. ( 3.10 ) together with the condition 
/? > uniquely determines (3. Then one finds 

lnP«-|ln(T/Ao), a ^ S[ (^13 + ^ - (3.11) 

One can verify that the convexity of 5*1 ensures the condi- 
tion a > if Ai < 1/t. Expression (3.11) determines the 
probability density fimction of ln(T/Ao). For the PDF 
of T we obtain 



a: 



a/2 



(3.12) 



Since A cx R^, the power-law distribution of A can be 
used to obtain the power-law distribution of the molecu- 
lar elongation, R [^. 

We see that the PDF is a power-law function with the 
exponent 1 -I- q;/2 that can be expressed via the entropy 
function Si . Since the precise form of Si is generally un- 
known, it is impossible to find the precise dependence of 
a on the parameters of the flow. However, some general 
properties can still be inferred. As Ai increases, i.e. when 
the Reynolds number increases, a decreases and tends to 
zero when Ai 1/t. One can easily establish the behav- 
ior of a in this region since th en t he approxim ation (|3.6|) 
is correct. Substituting Eq. into Eqs. (|3.10| , |3.11^ 

we obtain 



"= A 



T 



Note that the only characteristics of the flow entering this 
expression are the average value of pi and its dispersion. 

Power tail ( 3.12| ) means a slow decay of the proba- 
bility distribution of T = trA, which results in infi- 
nite values of its high moments. Namely, the moments 
(T") = / dTP(T)T" diverge if n > a/2. Moreover, the 
normalization integral / dT'P{T) converges at large T 



only if a > 0. Therefore an attempt to extend the pas- 
sive consideration to Ai > 1/t leads to the divergence 
of the normalization integral. It can be interpreted as 
the tendency of the polymer molecules to be stretched, 
i.e. the coil-stretch transition, and the breakdown of the 
passive approach. 

As we have seen, the n-th moment of the conforma- 
tion tensor A formally diverges at n > a/2. It signals 
the breakdown of the passive approach, i.e. the main 
contribution to the diverging moments comes from the 
configurations of the velocity such that the feedback of 
the polymers on the fiow cannot be disregarded. As ex- 
plained in the beginning of this section, the molecules 
can be considered as passive at 11 ^ i/Xi. An account 
of the back reaction of the polymers on the fiow leads 
to a much faster decrease of the PDF of 11 at 11 > i/Xi, 
which ensures that the moments have finite values (in the 
framework of the simple model presented in Appendix ^ 
one can find the precise form of the PDF). Let us estimate 
the value of the diverging moments taking the feedback 
into account. It is more convenient to discuss the mo- 
ments of n, w hich are proportional to t he moments of A 
(see Eq. {^)). If a < 2 then Eq. ( |3.12D modified by the 
cutoff at n ~ vXi gives 



(tr n) 



n; 



lyXi 



a-a/2 



(3.13) 



Note that (trll) Hq becau se we assumed that IIq <C 
vXi. The equations ( 2.1C , 3T3| ) show that the elastic con- 
tribution to the energy dissipation, T~"'^(trn), can be es- 
timated as T~^nQ^'^(i^Ai)^^"/^. It becomes comparable 
to the viscous contribution, ~ I'Xi, just at the point of 
the coil-stretch transition, where a — 0. 



C. Correlation Functions 

Here we investigate simultaneous many-point correla- 
tion functions of A. Let us start with the two-point cor- 
relation function, 



G2(r)-(T(to,ri+r)T(to,ri)) 



(3.14) 



where T = tr A. The value of A at a given point is deter- 
mined by the Lagrangian trajectory arriving at this point 
at t = to. Polymers separated by distances smaller than 
the viscous scale, ry, are stretched coherently, whereas at 
larger separations the correlation is largely lost. There- 
fore 7] is the correlat ion l ength of A. For distances r ^ rj 
the quantities T in ( 3.14 ) become weakly correlated and 
the correlation function tends to the product of aver- 
ages (T). Non-trivial correlations occur at the distances 
smaller than rj. The correlation function is a monotoni- 
cally decreasing function of the distance r. 

Note that G2 cannot be calculated in the framework 
of the passive approach if a < 2. Indeed, we can write 
G2 ~ (T)^ at r 77. Formal calculation of (T) using 



PDF gives an infinite result. The same is true for 
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G{0) = (T^). Since G2(r) is a decreasing function of r, 
it follows that G2 (r) is infinite in the framework of the 
passive approach. It means that the main contribution to 
the correlation function comes from configurations such 
that the back influence of the polymers on the flow is not 
small. On the other hand, if a > 4 the calculation of both 
(T)^ and (T^) in the passive framework gives a value of 
the order of Aq. Hence G2{r) ^ Aq for all r, which means 
that the main contribution is made by the configurations 
where the polymers are close to the equilibrium state. 

Let us consider the most interesting case 2 < a < 4. 
Then (T) ^ Aq is finite (and small), whereas (T^) is 



infinite if calculated using the "passive" PDF (B.12). It 
means that at small distances the back reaction must be 
taken into account, whereas at larger distances the pas- 
sive approach works well. If we calculate the correlation 
function in the framework of the passive approach, the 
result is valid for distances larger than some characteris- 
tic scale of the back reaction. 



To calculate corr elati on function (3.14) one can sub- 
stitute expression (2.17) for A{ri + r) and A{ri) and 
then average over ri which can be considered as aver- 
aging over space. The distance between Lagrangian tra- 
jectories terminating at the points ri + r and ri is a n 
increasing function of t. It can be found from Eq. ( [2.15 ): 



Sx{to - t) = Wito - t,to)r = W-^{to,to ^ t)r . (3.15) 



Expression ( 3.15| ) is correct provided \Sx\ < 77. Under the 
s ame condition the matrices W entering the expressions 
(2.17) for A(ri + r) and A(ri) are identical, as follows 
from Eq. (p. 14). Using the decomposition (3.1), rewrit- 



ten as W''^{tQ, to —t) — N'^ K^^ M'^ , and the inequalities 
gPi > qP2 ^ f,P3 obtain \5x\ « e^P^'r. 

At r <C 77 we have G2 ^ Aq, i.e. the main contribu- 
tion to G2 is due to the rare events when the product 
WW^ is anomalously large during a lo ng time. Then 
one can use the uniaxial approximation (3.7). The func- 
tions pi{t) in A(ri -f r) and A(ri) are identical as long as 
\5x\ < rj. When the separation 5x becomes larger than ry, 
the correlation between the Lagrangian trajectories be- 
comes weak. The contribution of this stage to G2 is given 
by the product of independent averages, (T)^ ~ Aq, and 
can be neglected. We conclude that the main contribu- 
tion to G comes from times when \5x\ < 77. Using the 
expression \5x\ ft! e~^^r we get 



dt exp(2pi - 2t/T) 



(3.16) 



where pz{s) = hi{r/r]). Since both pi and p^ ent er t he 
integral, to evaluate G2 one needs the joint PDF (3.4). 
To ensure a large value of G2, the function pi in Eq. 
.16) should first increase faster than t/r, and then re- 
turn to its average value Xit. Thus pi — t/r should have 
a maximum at < s. A vicinity of i = makes the 
main contribution to G2. In the absence of the con- 
straint t» < s the value of G2 grows exponentially as 



(1 



increases, which corresponds to a formally infinite value 
of the second moment at a < 4. Therefore the optimum 
is achieved at t» » s. It gives the estimate 



G2 - A^ J dpi exp(4/9i - 4s/r) 

'pi-Xis h-i{r/r]) - X3S 



X exp < —sS 



(3.17) 



where s is determined from the condition ^3(5) — ln(r/77). 
The integral can be calculated in the saddle-point ap- 
proximation with the saddle-point pi oc s. Next, one 
should optimize over s, which gives s oc ln(r/77). The 
proportionality coefficients depend on the form of S. The 
substitution of the optimal values gives 



G2 ^ Al [r^/r) 



i2 



(3.18) 



The exponent ^2 in Eq. ( p. 181) can be found if the precise 
form of the entropy function S is known. We observe 
that G2 ~ Aq if r ~ 77. It is natural since G2 at r ~ 77 
can be estimated as (T)^. We also see that G2 tends to 
infinity at r — *■ 0. This corresponds to a formally infinite 
value of the second moment. 

All the conclusions concerning the pair correlation 
function of T are valid for correlation functions of sepa- 
rate components of A t oo, which follows from the single- 
axis substitution (3.7). Indeed, n is uniformly dis- 
tributed over the unit sphere, reduces correlation func- 
tions of A to correlation functions of its trace T. 

Let us discuss the case a > 4. Then the main contri- 



bution to G2{r) (3.14) at small r is equal to the second 



moment (T^) ^ A^. One can examine the r-dependent 
correction to the second moment (T^) — G2(r). It can 
be done as above. The correction behaves as a positive 
power of r at r ^ 77. 

The proposed scheme can be generalized to higher- 
order correlation functions 



G„ = (A(ri)...A(r„)) 



(3.19) 



The behavior of G„ is similar to tha t of G2. If the mo- 
ment (T") calculated with the PDF (|3.12|) is infinite, the 
function G„ is a scaling function of the coordinates. The 
scaling exponent is negative, so the correlation function 
formally diverges at small distances. On the other hand, 
if the moment (T") is finite then the difference G„ — (A") 
scales with a positive exponent and is thus a small cor- 
rection to (A"). 

Since the moments of 11 are finite, we can assert that 
the growth of the correlation functions of 11 observed at 
fusing points in Eq. (3.19) has to be saturated. F or ex - 
ample, at 2 < a < 4 the pair correlation function (3.14) 



saturates at G2 ~ ^0-^0 ' 



A2n[;/'-'(ivAi)2-"/2. 



One can say that the back reaction regularizes the corre- 
lation functions at small scales. 
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IV. STRONG BACK REACTION 



Here we consider the dynamics of the polymer solu- 
tions above the coil-stretch transition, when the Reynolds 
number exceeds a critical value. Depending on the con- 
centration of polymer molecules two situations are pos- 
sible. If the concentration is small, the elastic stresses 
are small in comparison with the viscous stresses. Then 
the polymers are stretched to their maximal elonga- 
tion, i?niaxj and the properties of the fluid do not dif- 
fer significantly from those of the pure solvent. Below 
we consider the more interesting case where the con- 
centration of polymers is large enough, so that elastic 
stresses can be larger than the viscous stresses. Then 
the feedback due to the polymers substantially modifies 
the flow. The condition for the existence of the back 
reaction regime is Xlmax S> v jr^ where Ilmax is the maxi- 
mal value of the elastic stress tensor. It can be expressed 
in terms of microscopic parameters and the concentra- 
tion a^he polymers as Ilmax = Ko''T'Q~^ ^max (see sub- 
sect 



the value of the gradient cannot exceed l/r, the velocity 



,11 A). Using estimates for the microscopic parame- 
ters proposed in Ref. ]T3[ | one can rewrite the condition 

as n> (i?o^max)"^- 

Whereas in pure solvent typical gradients of the ve- 
locity grow unlimited as the Reynolds number increases, 
in polymer solutions the balance of inertial and elastic 
degrees of freedom fixes the characteristic value of the 
gradient at l/r. Indeed, if the instantaneous velocity 
gradient exceeds 1/t, it extends the polymers, so that 
the elastic stress grows and damps the gradient. On the 
other hand, if the velocity gradient is much smaller than 
1/t, the molecules contract and do not produce effect on 
the flow. Then the velocity gradients tend to grow to the 
value characteristic of the pure solvent, which is larger 
than l/r above the transition. This explains the steady 
state realized above the transition. We now establish 
some general properties of this steady state. 

Turbulence of Newtonian fluids can be characterized 
by two length scales: the integral scale, L, and the dissi- 
pation scale, 77. Energy pumped at the integral scale cas- 
cades without dissipation from larger to smaller eddies 
(coherent motions of the fluid) in the range rj < r < L, 
called the inertial interval. Velocity difference between 
two points separated by the distance r from the inertial 
range diminishes slower than r, so that the characteristic 
value of the velocity gradient at the scale r grows down- 
scales reaching a maximum at r ~ 77 |2^] . If Vt /L 1, 
where V is the velocity at the integral scale, the gradient 
related to large eddies is smaller than r^^. Therefore, 
large eddies do not excite polymers, which means that 
the elastic stress tensor is not correlated at these scales. 
Since only coherent excitations of the elastic stress tensor 
can influence the velocity, we conclude that the elasticity 
is negligible for large eddies. The interaction of inertial 
and elastic degrees of freedom becomes essential at the 
scale r^t, determined from the condition Vw ~ l/r. The 
fluctuations of 11 are correlated over the scale . Because 



difference scales linearly with r at r 



< 



i.e. the flow 



is smooth. Near the coil-strctch transition characteristic 
velocity gradient is determined by the viscous scale and 
is of the order of l/r, hence 77. As the Reynolds 

number increases, velocity fluctuations increase, so that 
the scale r, grows. Thus above the coil-stretch transition 
a new scale, r*, separating the inertial and viscoelastic 
intervals arises. It is of the order of rj near the transition 
and grows as the energy input increases. 

N ear the transition the viscous and elastic terms in Eq. 
( ^.8[ ) are of the same order, which gives 11 ~ v/t. For 
dilute solutions v/t is much larger than Hq, there fore all 
the terms in the energy balance equation (2.1C) are of 
the order of v/t'^ near the transition. The energy pump- 
ing rate per unit mass, e, can be estimated as V^/L. 
Equating it to the dissipation rate, estimated as v/t'^, 
one finds the estimate Rcc = [L^/(z/r)]^/'^ for the value 
of the Reynolds number at the transition. As the energy 
input increases the energy dissipation rate due to viscos- 
ity, 7^(Vv)^, remains of the order of i^/r^. Therefore far 
above the transition the principal part of the energy is 
dissipated by the polymer relaxation. Then the viscous 
term in Eq. ( 2.1C| ) can be neglected and we obtain 



(tr n) 



(4.1) 



We conclude that the energy is d issipated mainly by the 
elastic relaxation. Relation (4.1) means that the typical 



value of n grows as the energy input increases, which can 
be interpreted as the increase in the effective ("elastic") 
viscosity. It is defined as the proportionality coefficient 
between the polymer stress tensor 11 and the strain ten- 
sor ViVk + Vfetij, which remains of the order l/r. Using 
Eq. ( 4T ) one can estimate the ratio of th e elastic term to 
the non- linear inertial term in Eq. (2.8) as Vt/L <C 1, 
which shows that the elasticity is indeed negligible at 
large scales. 

The strong interaction between the elastic and iner- 
tial degrees of freedom imposes a restriction on the La- 
grangian statistics of velocity. To demonstrate it, ob- 
serve that Eq. (2_^) gets simplified under the condition 
n ^ Ho satisfied in the strong back reaction regime. 
Neglect ing the terms proportional to Aq and IIq in Eqs. 
( p^ 2.7) we obtain 



dtn,j + (t;V)% = nfcj Vfcw, -f- u.kVkVj - -u,j , (4.2) 



Expressing the solution of Eq. (4.2) in terms of the 
Lagrangia n quantities x and W introduced by Eqs. 
(|ti|,|t|) we obtain 



n(i, r) = W{t, 0, r)n[0, x{0, r)]W^{t, 0, r)e 



-2t/T 



(4.3) 



The Lagrangian correlation time at the scale is r. 
Therefore at i ^ r the eigen-values of W are strongly 
separated so that 11 is uniaxial: 



Ily = rLiUj trll , 



(4.4) 
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where n is a unit vector. Then Eq. (^^) gives 
2t 

2pi = — +ln[trn(t,r)]-ln[trn(0,a;(0,r))]. (4.5) 

T 

The stationarity of 11 impHes that pi — t/r has a sta- 
tionary distribution. In particular, we conclude that the 
principal Lyapunov exponent Ai of the flow is equal to 
1/r exactly. The stationarity of pi — t/r is very different 
from the situation for the Newtonian fluids, described by 
Eq. (3^). The reason is the anticorrelations in the tem- 
poral dynamics of a due to its interaction with 11 which 
were qualitatively described in the beginning of the sec- 
tion. They lead to vanishing dispersion A of pi — t/r, 
A = / di((CTii(t)(Tii(0))), which is non-zero for Newto- 
nian turbulence. 

Averaging Eq. (|4.2|) one can obtain the exact relation 



(n,;fc(r)Vfcz;,(r)) = 



(trn) 



(4.6) 



Consider now (nifc(r)V/cWi(r')) as a function of the sepa- 
ration I = \i — r'\. Its value at Z = is given by Eq. (4.(;) 



and can be shown to be much larger than the value at the 
pumping scale, I ^ L. Indeed, consider the correlation 
function averaged over a ball of size L centered at r, i.e. 
ljlik{r) Jy dr'WkVi{r'))/V. The velocity gradient aver- 
aged over th e sc ale L is determined by the external forces. 
Using Eq. ( 11 ) one can estimate the value of the aver- 
aged correlation function as cVt/L ^ e. It follows that 
(nifc(r)V/jWi(r')) decays at scales larger than r,. Below 
the fluctuations of 11 and Vu are strongly correlated. 
The decay of the correlation function at < r < L can 
be used to derive the Kolmogorov's four fifths law |2^] 
at these scales. The latter states that the third order 
longitudinal structure function is equal to — 4e^/5 in the 
inertial interval. All the above conclusions are in agree- 
ment with the general picture presented in the beginning 
of the section. 



Expression (4T) gives the typical value of the stress 
tensor. As we argued above, the fluctuations with 11 ^ 
er relax rapidly due to the back reaction, which leads 
to a fast decrease of the PDF of 11 at 11 ^ er. On the 
other hand, the probability to have 11 ^ er is also small. 
The rough details of the behavior of the PDF can be 
understood on the basis of a simple model presented in 
Appendix The solution of the model shows that the 
PDF of n has an exponential tail at large values of 11 
and power-law behavior at small values of 11. We be- 
lieve that a similar qualitativ e behav ior is realized for 
the stress described by Eqs. (2^,12). The model also 



explicitly demonstrates the finite value of the Lagrangian 
correlation time of 11 and Vf . This property holds de- 
spite a strong modification of the Lagrangian dynamics 
due to the back reaction. 

Note, that the concentration of the pol yme r molecules, 
n, does not enter system of equations (2.8,4.2). Therefore 
the dynamics of polymer solutions with different values 
of n will be identical in the strong back reaction regime. 



Moreover, using the equation dtu + (t;V)n = for the 
concentration, it is possible to show that Eq. (12) is 
also valid for inhomogeneous in space n. Thus the hydro- 
dynamic properties of spatially inhomogeneous solutions 
do not differ from the homogeneous ones. This assertion 
holds if local n is large enough for nniax(?^) to be larger 
than the local value of 11 prescribed by the dynamics. 



The uniaxia l form (4.4) of the tensor 11 allows one to 
rewrite Eqs. (2.8 4.2) in the form similar to the equa- 
tions of the magnetic hydrodynamics. The field n Vtrll 
satisfies the induction equation with linear damping. In 
addition one can show that the field is solenoidal. This 
analogy helps understand the dynamics of fluctuations 
at the scales r <C r*, which occur on the background 
of the relatively slow stresses excited at r ~ r* . These 
small-scale fluctuations are elastic waves similar to the 
Alfven waves propagating in the presence of a large-scale 
magnetic fleld in magnetic hydrodynamics |26 2^. The 
dispersion relation of the waves is w = ky/trll. Thus the 
velocity of these waves is -y/trTT which can be estimated 
as \/eT. There exist two mechanisms of the elastic waves 
attenuation: polymer relaxation and viscous dissipation. 
The flrst mechanism leads to the scale-independent at- 
tenuation , which is smaller than the frequency, at 
fcr» 3> 1. The second mechanism produces the attenua- 
tion ~ vk^ which is much smaller than the frequency for 
fcyy* 1 where 77* = z/(eT)~^/^. Thus the elastic waves 



are well-deflncd in the interval r^^ ^ k ^ rj^^ . 

Our equations are valid as long as 11 <C Ilmax- The 
relation (4.1) allows us to reformulate this condition as 
e <C ninax/T. Another limitation of our scheme is related 
to the inequality i? <C , under which the flow is smooth 



at the scale R. Using Eqs. (|2.5|), (p7|), and (4.1) one can 
write the estimate ^ QeT{K{)n}^ for the typical size 
of a polymer molecule, R. 

Let us estimate the parameters introduced above 
within the framework of Kolmogorov's theory (K41) 
psf . Though the theory is, rigorously, incorrect [ p5| , 
it is satisfactory for rough estimates. The characteris- 
tic velocity difference, 5rV, between two points separated 
by the distance r from the inertial interval is given in 
K41 by (er)^/'^, where e is the energy input. Writing 



|Vi;| 



,l/3„-2/3 



finds 



In the 



K41 theory the condition R can be rewritten as 

g(KQnT'^)~^ <C 1. Note that in the framework of K41 
theory the ratio R/r^ is independent of the Reynolds 
number. 

Our analysis assumes that the characteristic size of the 
molecules, R, is much smaller than their maximal size, 
-Rmax- As Re increases, the typical elongations eventu- 
ally become of the order of i?max, and further elongation 
becomes impossible. In this case the molecules behave as 
hard rods, modifying the effective viscosity of the fluid 
p9[ . Therefore at large enough Re we return to the case 
of Newtonian fluid. However, this regime is expected to 
be unstable because polymer molecules are intensively 
destroyed by strong flows. 
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We have shown that in the steady state the velocity 
gradients in the bulk do not exceed t~^. Consider now 
the situation where the boundary forces tend to produce 
gradients larger than 1/t at r ~ L. Then the elastic reac- 
tion should lead to formation of a boundary layer where 
the value of velocity gradient diminishes from the value 
imposed by the forcing to the value r^^ in the bulk. Then 
r* ~ L, i.e. the inertial range and energy cascade are ab- 
sent. This situation is similar to the elastic turbulence 
regime 

Finally, let us consider the role of other modes of the 
polymer molecules relaxation. They are characterized 
by the relaxation times Ti < r. We have shown that 
the interaction of the fluid with the principal relaxation 
mode fixes the value of the principal Lyapunov expo- 
nent at Ai = T~^. The inequality Ait^ < 1 then implies 
that other modes are always only weakly excited by the 
flow, so the interaction is fully determined by the soft- 
est relaxation mode. We conclude that the equations 



( 4^ ) based on the single relaxation mode approxima- 
tion correctly describe the solution hydrodynamics above 
the coil-stretch transition. 



V. CONCLUSION 

We have examined properties of turbulence in dilute 
polymer solutions. Our results support the theory of 
Lumley ||] who argued the existence of the coil-stretch 
transition in turbulent flows, which occurs at a criti- 
cal Reynolds number, Rcc. The polymer molecules are 
typically weakly stretched, so that their elasticity only 
weakly influences the flow in the regime realized below 
the transition, at Re < Rcc. At Re > Rec the polymer 
molecules are substantially stretched and strongly mod- 
ify the small-scale flow. 

At Re < Rcc the polymer molecules are passively ad- 
vected and stretched by the flow. This regime occurs 
under the condition Air < 1 where r is the polymer 
relaxation time and Ai is the principal Lyapunov expo- 
nent. The Lyapunov exponent is defined as the logarith- 
mic rate of the divergence of nearby Lagrangian trajec- 
tories and can be estimated as the inverse turnover time 
at the viscous scale of turbulence. The majority of the 
molecules in this regime fluctuates near the equilibrium. 
There also exists a small number of strongly elongated 
molecules which appear due to rare large fluctuations in 
the rate of strain. Even though the number of substan- 
tially elongated molecules is small, they may be relevant 
in some situations due to the relatively slow power-law 
decrease of th e pro bability density function of elongations 
of molecules ( 3.12 ). 

In the second regime, at Re > Rcc, most of the 
molecules are substantially elongated. It leads to a strong 
interaction between the elasticity and the flow, which 
modifies the flow below the scale . At r > r* the prop- 
erties of turbulence are the same as in Newtonian fluids. 



The energy cascades downscales from the pumping scale 
and dissipates due to polymer relaxation at r ~ . The 
scale can be considered as a new dissipation scale. The 
flow is smooth at r r* with the Lyapunov exponent Ai 
fixed at the value 1/r by the interaction. 

The smoothness of the flow at r < leads to the con- 
clusion that the velocity spectrum E{k) decreases faster 
than at fcr* ^ 1. The precise form of E{k) in this in- 
terval is related to the elastic waves propagating at these 
scales. As both spectral transfer time and the decay time 
are scale-independent one can expect a power-law spec- 
trum. 

The properties of the polymer statistics near Rec were 
examined numerically b y Kronj ager and Eckhardt ISOfl in 



the framework of Eqs. (2.6,2.8). The results indicate the 



power PDF tail for the polymer elongations at Re < RCc 
and a substantial modification of the PDF at Re > Rec, 
in agreement with our results. 

Let us discuss implications of our results for the drag 
reduction. A description of the experimental situation 
can be found in the works It has been observed 

that the onset of the drag reduction at increasing Re 
depends on the concentration of the polymer molecules, 
whereas asymptotically the friction force falls on a curve 
which is independent of the concentration. This curve is 
usually referred to as MDR (maximum drag reduction) 
asymptote. A discussion of the MDR can be found in 
the recent work [^ . A natural explanation of the n- 
independence of the MDR asympt ote can be formulated 
in the framework of Eqs. (^.^,4.2) describing the strong 
back reaction regime. Indeed, the system contains no n- 
dependent parameters. The n-dependence of the onset 
can also be explained in our scheme. The drag force is 
formed in the boundary layer which has a complicated 
structure Q. Whereas gradients of the average veloc- 
ity grow toward a wall, the amplitude of the velocity 
fluctuations decreases. Therefore one can expect that 
the polymer molecules are strongly extended in bulk and 
weakly extended near the walls. Then the structure of 
the boundary layer will be sensitive to the polymer con- 
centration n. The situation corresponds to the transient 
regime (which is sensitive to the polymer concentration) 
from the Newtonian behavior to the MDR asymptote. 
The asymptote itself corresponds to the case when the 
polymer molecules are strongly extended everywhere. 

A striking property of polymer solutions is the so-called 
elastic turbulence regime discovered by Groisman and 
Steinberg jlj. It is a chaotic state which is realized at 
small Reynolds numbers Re. Its existence is made possi- 
ble by the large value of Weissenberg number Wi — tV/ L 
which implies a strong non-linearity of the system. This 
state can also be investigated in the framework of our 
scheme. The results will be published elsewhere. 

Let us give numerical values of parameters appearing 
in our theory for a typical experimental arrangement. For 
the number of monomers 10^ — 10'' one has Rq ~ 10~^ 
cm, i?inax ^ 10^^ cm, and r ~ 10~^ — 10~^ s. Then using 
ric ^ (R^^^Rq)^^ one can obtain 0.1 ppm for the concen- 
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tration, ric, below which polymers have no effect on the 
flow. Let us assume that the polymer concentration is 10 
ppm, the integral length is L '--^ 10 cm and take the water 
viscosity, i/ ~ 10~^ cw? j s. Then the critical Reynolds 
number Rcc ~ {L"^ / {vt)Y^'^ is of the order of 10^. 
Above the coil-stretch transition the characteristic si ze 
of polymers is given hy R ^ Rq ^Re^ gri^^ / {kBTnL'^) ^ 

10"^i?oR.e^''^- We obtain that R - lO^i^o in the vicinity 
of the transition, which is in agreement with the assump- 
tion Rq <^ R Rmax- Using the Kolmogorov's estimate 
77 ^ LRe"^^'* wc find that at the transition 77 ^ 10~^ 
cm, which is of the order of i?max- Thus the assump- 
tion i? ^ 77 is satisfied. These estimates seem to fit the 
existing experimental data. 
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APPENDIX A: 

Here we consider the statistical properties of the inte- 
gral 



I 



dt exp 



dt'^{t') 



(Al) 



where ^(t) is a random process with a finite correlation 
time and a negative average ^0 < 0- For the purpose 
it is convenient to introduce the auxiliary object 



lit) 



Due to the stationarity of ^ the statistics of I{t), it is 
independent of t. Separating the integration interval one 
finds the relation 







/ df exp 









I{t - 5t) = I{t) exp 



t-St 



^{t')dt' 







I dt' exp 


f mdt" 


h-st 


Jt-St 



It follows that 



Inl{t-St) = liil{t)+ [ ^{t')dt' 

Jt-St 



hW 1 



It-st dt' exp 


'it'-st ^indt"' 


I{t) exp 


'listmdt' 





(A2) 



If I{t) is sufficiently large, it is possible to neglect the 
last term on the left-hand side. The exact condition is 
formulated below. Observe that ln/(t) depends on the 
values of the noise at times larger than t so that the sec- 
ond term is independent of the first provided 6t t^. 
Therefore the probability distribution function, P{z), of 
z{t — St) = \nl{t — St) is given by the convolution of the 
distributions of z{t) (which is also equal to P{z)) and 
J^^g^^{t')dt'. The latter has a probability function simi- 
lar to the one (|3.5|). We thus obtain the integral equation 



P(z) 



dz' 



%/27rfA 



P(z') exp 



'StSf 



^o5t 



St 



where ^0 = (0 ^-i^d is the entropy function charac- 
terizing S^{t). Since the kernel of the integral operator 
depends on the difference z — z' only, the solution of this 
equation is P{z) cx exp[— az]. We obtain the following 
expression for the tail of the PDF of / = 



P(/) 



Here a is determined from the condition 



dx 



V27rfA 



exp 



StS. 



X - ^pSt 
St 



= 1. 



(A3) 



(A4) 



The solution a = should be rejected. This integral can 
be evaluated by the saddle-point method so that its value 
is determined by the maximum of the exponent. Taking 
its value a t x — ^oSt we conclude that in order to satisfy 
condition (A4), a and ^0 niust have different signs. Thus 
if ^0 < 0, then a > and the normalization integral for 
PDF (A3) converges at / — > 00. On the other hand, if 
^0 > 0, there exists no well-defined distribution of /. 

The two equations that implicitly define a are given 
by the saddle-point condition 

a = S'^{(3 - ^0), 

where /3 is the saddle-point value of {z — z')/St, and the 
condition 



which follows from the condition that the integral (A4) is 
equal to 1. One should reject the formal solution /? = <^o 
of these equations corresponding to a = 0. It is easy 
to see that (3 is positive together with a. Now we may 
formulate the condition for the applicability of the power 
tail. It is valid provided the third term in Eq. ( |A^ ) is 
indeed much smaller than the second for those I{t) that 
determine the PDF of I{t — St). From z~ z' = (3St it fol- 
lows that I{t) = I{t — St) exp[—(3St], so that we arrive at 
the condition / exp[f}St]/{P'^St) (^ is estimated as f3). 
The increment St constrained by the condition St ^ t^. 
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There are two cases to be considered. If ^ one can 
use the c hoic e minimizing the above ratio 6t ~ so 
that Eq. ( |A^ ) is vahd for I ^ f3~^. In the opposite case 
f3~^ < the power tail is vahd for / ^ exp[/3T^]/(/3^T^). 

At smah ^0 one can use quadratic expansion for 
S^{x) ~ x^/{2A) which gives 



P = -£.0, Oi^-—£,o. 



(A5) 



The entropy function becomes quadratic in the hmit 

^ 0. Thus the expression (A5) is vahd for any 
in the case of a short-correlated process ^. 



APPENDIX B: MODEL OF THE BACK 
REACTION 

Let us introduce a simple model that captures the most 
robust features of the interaction between elastic and in- 
ertial degrees of freedom. The model is formulated in 
terms of the system of equations for two variables a and 
X. The equations are 



dx 

— ^ ax + xq 
at 



(7 = -x + ^, (Bl) 

mtinh))) ^ 26{h ~ , (B2) 



where double brackets denote the irreducible part of the 
correlation function. The variable a models the rate-of- 
strain subtracted by 1/t and x models the elastic stress. 
The time derivative in Eq. (Bl) represents the full deriva- 
tive dt + I'V, i.e. we consider the Lagrangian dynam- 
ics. The product ax represents the combination on the 
right-hand side of Eq. (4.2) and xq stands for Aq in Eq. 
(^■6\) . The second equation in system (Bl) represents 
the Navier-Stokes equation (2^). Since we consider dy- 
namics at the scale r*, all the spatial derivatives can be 
estimated as 1/r*. The term —x describes the back reac- 
tion and ^ models the influence of larger scales, exciting 
the motion at r '-^ r*. The average a = {£) is negative 
below the coil-stretch transition and positive above. If 
a > the t erm xq on the right-hand side of the first 
equation in (Bl) can be disregarded. 



Starting from the system of equations (B1,B2) one can 
derive the Fokker-Planck equation for the PDF of x: 



dtV = [xd^ (xV)] - d., [{xa + (a - x)x)V] 



The normalized stationary solution of Eq. (B3) is 



V^{x) 



exp 



(B3) 



(B4) 



is the same value as without the back reaction. However, 
the back reaction is important for rare events when ^ is 
large. The interaction leads to the exponentially decay- 
ing tail which makes all the moments of x finite. This 
corresponds to the picture presented in the main body of 
the text. Note that the power tail is universal, i.e. force- 
independent, whereas the exponential tail is an artifact 
of a zero correlation time of ^ . 

In the case a > 0, i.e. above the transition, the limit 
Xq — !■ is regular. One obtains 



Vo{x) 



exp(- 



T{a) 



(B5) 



where r(a:) is the Euler F-function. We observe that all 
the positive moments of x exist, because the back reac- 
tion stops the growth of x. The average value of x is 
given by (x) = a so that {a) = 0. These facts correspond 
to the statements (11) — er and Ai = from Sec. IV 



Let us now investigate non-simultaneous correlation 
functions of x above the transition, i.e. when a > 0. 
Then we can assume xo = 0. We need the Green func- 
tion Q{t, X, y) of Eq. (B5), which satisfies 



dtQ-dx [xdx [xQ)] + dx [{a - x)xQ]^5{t)5{x~y) , (B6) 

with the condition Q{t < Q) = 0. The Green func- 
tion should be regular at x = and decrease faster 
than any power of x at x ^ oo. Using Q one can find 
non-simultaneous correlation functions of x in the steady 
state: 

(/i[x(t)]/2[x(0)]) = j dxdyVoiy)g{t,x,y)hix)f2iy), 

where /i and /2 are arbitrary functions and Vo{x) is de- 
fined by Eq. ( p35| ). 

The Laplace transform of g{t, x, y) is 

poo 

G{X,x,y)= dt exp(~Xt)g{t,x,y). (B7) 

"'0 

It satisfies the equation 

AG = dx [xdx (xG)] - dx [(a - x)xG] + 5{x - y) , (B8) 



following from Eqs. (B6,B7). The solution of Eq. ( pjq ) 
can be expressed in terms of two independent solutions 
of the homogeneous equation 



AGi,2 = dx [xdx ixGi,2)] ~ dx [(a - x)xGi 



where Z = 2x'^^^ Ka{2y/^) is the normahzation factor. Two independent solutions of Eq. (§9|) are 



Here Ka is the MacDonald function. 

At a < 0, which corresponds to system (2.6-p^) below 
the transition, the properties of a are only insignificantly 
modified by the interaction with the variable x. For ex- 
ample, in the limit xq <C 1 one finds that (a) = a, which 



Gi = x''^e-''F{-k2 - 1, fci - /C2 + 1, x) , 
Ga = x'''^e"''^'(-fci - 1,A:2 - fci + l,x) , 



(B9) 



(BIO) 
(Bll) 



where 
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r(i-/3) 



r(Q; - p + 



r(/3-i) 
r(a) ' 



:i"''F(a-/3 + l,2-/3,x). 



(B12) 



and F{a,(3,x) is the confluent hypergeometric function 
[ p3[ . The functions Gi and 6*2 satisfy the boundary con- 
ditions at X = and a; ^ oo correspondingly. The pa- 
rameters ki 2 are 



2 ± Va2 -f 4A 



(B13) 



Matching functions (B10,B11) at x = y one obtains the 
Green's function 



n-1 



G{X,x,y)= , „ , j^f, ,x 
Va2 -I- 4A i^{ki - k2) 

+9{x ~ y)Gi{y)G2{x)] , 



[9{y ~ x)Gi{x)G2{y) 
(B14) 



where 9{x) is t he H eaviside step function. 

Expression (B14) allows one to establish analytical 
properties of G as a function of A. The function is ana- 
lytic in the half-plane ReA > 0. There is a branch point 
at A = —0^/4 with the cut going along the axis ImA = 
from the branch point to — oo. In addition to the branch 
point there also exist poles located at A = n{n — a) 
where n is an integer number (including zero), such that 
n < a/2. Thus the poles lie between the origin and the 
branch point. One can easily find the pole contribution 
to G near A = corresponding to n = 0. Using Eq. 



(B14) we obtain 



G{X,x,y) = 



exp(— a;) 
r(a)A 



(B15) 



where dots mean terms regular in A. 

The Green function Q is expressed via its Laplace 
transform as 



Q{t,x,y) 



A+ioo 

— cxp{\t)G{\,x,y), (B16) 

A—ioc 



27ri 



where the integration contour lies on the right of all sin- 
gularities of G(A). Shifting the integration contour in 
Eq. (B16) to the left w e fin d that the first term in the 
right-hand si de of Eq. ( B18 ) is reproduced by the pole 
contribution (B15) and then 



— eM^t)G{X,x,y), (B17) 

where e > is a small number and Qi is defined by the 
relation 



Gi{t,x,y) = gi{t,x,y) - Vo{x) . 



(B18) 



Shifting the integration contour in Eq. (B17) to the left 
we encounter the branch point A = —ayAii a < 2 or 



the pole A = l — aifa>2. Therefore at large t the t- 
dependence of Qi is cxp{—a^t/4:) if a < 2 and exp[(l — a)<:] 
if a > 2. Since 



{fi[x{t)]Mx{Q)]) = if,) ih) 

\dxdyVQ{y)Gi{t,x,y)fi{x)f2{y) ■ 



(B19) 



the correlations of x decay exponentially in time. 

One can write the asymptote of Qi {t, x, y) at large t ex- 
plicitly. Let us first c onsider a < 2. Then the integration 
contour in Eq. (B17) can be deformed into a curve going 
around the cut starting from A — —0^/4. Calculating the 
jump on the cut we get the expression 



Giit,x,y) 



2 ln(?;/x)ci + c 



16i3/2, 



1 a/2-1 -a/2 - 
X ' y ' e 



(B20) 



valid for small x, y. Here ci = 2ilj{l) — ip{—a/2) and 'ip{z) 
is the logarithmic derivative of the F-function. At a > 2 
with the exponential accuracy the function Qi(t,x,y) is 
given by the residue at A = 1 — a: 



Giit,x,y) 



{a-l)F{-l,a-l,x)F{-l,a-l,y) 



T{a-2)y 



Finite correlation time of a follows from {a{t)a{0)) = 
—df(\nx{t)h\x{0)). The expectation value of x{t) = 
exp[Jp (T{t')dt'] grows exponentially with time at times 
much larger than the correlation time of a unless 
j dt{a{t)a{Q)) = 0. Thus the back reaction stops the 
growth of x and gives rise to this peculiarity of the statis- 
tics of a. 
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